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Abstract. Given a vector space of microscopic quantum observables, density 
functional theory is formulated on its dual space. A generalized Hohenberg-Kohn 
theorem and the existence of the universal energy functional in the dual space are 
proven. In this context ordinary density functional theory corresponds to the space 
of one-body multiplication operators. When the operators close under commutation 
to form a Lie algebra, the energy functional defines a Hamiltonian dynamical system 
on the coadjoint orbits in the algebra's dual space. The enhanced density functional 
theory provides a new method for deriving the group theoretic Hamiltonian on the 
coadjoint orbits from the exact microscopic Hamiltonian. 
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1. INTRODUCTION 

Density functional theory, in its simplest form, is a theory for nondegenerate ground 
states that replaces the many particle antisymmetrized wave function by its density p{r). 
The Hohenberg-Kohn theorem shows that there exists a universal energy functional 
minimized by the density of the exact ground state wave function Jl], @]. In addition 
to ubiquitous applications in quantum chemistry this many body theory solves 
problems in solid state |§, H and nuclear physics 0. 

This letter's aim is to generalize density functional theory to a framework where the 
densities are elements of the dual space to a vector space of microscopic observables. The 
Hohenberg-Kohn theorem is proven in this context, and the existence of the universal 
energy functional on the dual space is established. If the vector space is also a Lie 
algebra, then a group theoretic model is defined naturally on the coadjoint orbits 
of the algebra's dual space. The generalized Hohenberg-Kohn theorem enables the 
derivation of the group theoretic Hamiltonian from the microscopic Hamiltonian on 
Fock space. The problem of constructing an algebraic model Hamiltonian from the 
microscopic interaction is an old one in group theory with many potential applications 
to the description of complex many-body systems. 

Let TC denote the Hilbert space of iV-fermion antisymmetrized state vectors, V{H) 
its complex projective space, and H the self-adjoint Hamiltonian operator on Ti. 
Suppose g is a real linear vector space of hermitian operators defined on H, and g* 
is its dual space, i.e., the vector space of real- valued linear functional on g. Denote the 
pairing between a dual element p E q* and an operator X G g by (p, X) = p{X). Define 
the moment map M : V{T-C) — * g* from the projective space of states to the dual space 
by 

(M(,),..).fM. 

where G 7i, X G g, M(0) G g*. A brief history of the moment map is given on p. 327 
of [0. Define the energy mapping F : P(7i) R from the projective space to the real 
numbers by 

and for each fixed X G g, define Fx : V{T-C) — R by 

Fxi<P) :=F(0) + (M(0),X). (3) 

For most physical systems only a few quantum states are observed or are of direct 
interest, e.g., the ground state, rotational and vibrational bands, collective states, and 
resonances. Although the microscopic theory on the Hilbert space Ti is complete, it may 
obfuscate the true physics because all degrees of freedom are treated democratically. 
By selecting out the most relevant degrees of freedom, a discriminating model based 
on the observables in g can elucidate the physics and simplify the mathematics. The 
moment mapping is an indispensable tool for creating a tractable theoretical description 
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of microscopic many-body systems. The density p = M{(f)) retains only part of the entire 
information about the system that the wave function (p carries, but a very important 
part - the expectations of the observables that span g. 

When g is a subspace of real- valued measurable functions v{r) on R^, and g* is the 
space of functions 



0* = p : ^ R 



\p{r)\d^r < oo L (4) 



conventional density functional theory is recovered. The elements of g are regarded 
as hermitian multiplication operators on 7i by associating the one-body operator 
J2n=i^i^n) with the function v{r) E g. The pairing is given by {p,v) = f p{r)v{r)d^r. 
If (j) is normalized, then the moment map M{(j)) = p is 

p(r) = N [ \(p{r,f2,...,fN)\'^d^r2...d^rN. (5) 



The functional F^{(f)) is the expectation of the energy including the one-body external 
potential Sn'^l^n)- 

Von Neumann's density matrix formulation of quantum mechanics is included in 
the formalism presented here by letting g be the space of all (not just one-body) 
hermitian operators on Ti. The dual space g* may be identified with a space of hermitian 
operators where the pairing is defined by {p,X) = tr{pX). The moment map is given 
by M{(j)) = I ) ( I if G 7^ is normalized. In this case, the moment map is one-to-one, 
and the microscopic theory is physically and mathematically equivalent to the density 
matrix expression of quantum mechanics. 

Another example is provided by cranking Hamiltonians or Routhians for which g 
and g* are isomorphic to the real three-dimensional vector space R^. The hermitian 
operator corresponding to cJ G q is u ■ J, where J denotes the one-body angular 
momentum operators on Ti. The pairing is the dot product, {p,uj) = p- u, and the 
moment map is M{(f)) = p, where 

= w 

Note that Fi^{(j)) is the expectation of the cranking Hamiltonian. 

In the case of the cranked Nilsson model, the space g consists of all hermitian 
operators formed from linear combinations of the angular momentum J and the 
quadrupole operator Q^^^ = J2n ''^n^/i'^K^n) ■ This eight dimensional space g is the 
rotational model Lie algebra rot{3). The full Nilsson model expands g to a 10 
dimensional space that includes the spin-orbit interaction and the square of the orbital 
angular momentum. Several more examples are produced by taking the vector space of 
operators to be the Elliott su{3) algebra the Ginocchio so(8) fermion algebra 
the general collective model algebra gcm{3) the symplectic model algebra sp{3, R) 
|TT|, and the unitary algebra of all one-body operators. In many of these examples the 
vector space g is a Lie algebra with respect to operator commutation. In §3 it is shown 
that when this additional algebraic structure is present, a mean field theory is naturally 
defined on the coadjoint orbit surfaces in the dual space g* . 
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2. HOHENBERG-KOHN-LEVY 

In this section the fundamental Hohenberg-Kohn theorem is proven in the extended 
framework of an arbitrary vector space of operators g on Fock space. The Levy 
constrained-search formulation is also established in this general case. 

Define the set of states S consisting of every G V{T-C) that is an absolute 
nondegenerate minimum of Fx for some X G 5. Set C equal to the image of S under 
the moment map, C = M{S) C g*. 

Theorem 1. The moment map M : S ^ C is an injection. 

Proof: The proof of this theorem is similar to the original Hohenberg-Kohn argument 
imi. Suppose, to the contrary, that there are two distinct states 0i,02 G 5, 0i 7^ 02, 
yet M(0i) = M(02)- There are operators Xi,X2 G g such that 0i is an absolute 
nondegenerate minimum for Fx-^ and 02 is an absolute nondegenerate minimum for 
Fx2- Since 0i 7^ 02, the nondegeneracy of the minima implies 

E, = Fx,{4>,) < Fx,{^2) 
E2 = FxM<FxM. 

g is a vector space, so Xi — X2 G g and 

El < Fx, (02) = F(02) + (M(02), (Xi - X2) + X2) 

= i^X,(02) + (M(02),(Xi-X2)) 

or Ei<E2 + (M(02), (Xi - X2)). 
By repeating the argument for interchanged 1 and 2, 

E2<El + (M(0i),(X2-Xi)). 

Adding the two equations, and recalling that M(0i) = M(02), leads to the contradiction 

Ei + E2<E2 + El. 

Hence, for two distinct states in S the moment map yields two distinct densities in C. 

The theorem depends only on the vector space structure of g, not on any possible 
additional algebraic structure. It is a corollary that for each X G g, the mapping 
Fx = Fx o : C R is well-defined. The state minimizes Fx on S if and only if 
p = M{(f)) minimizes Fx on C. 

A density p is said to be v- or, in the general setting, X-representable if p G C = 
M{S). The energy functional must be minimized on the set of X-representable densities. 
If the X-representable densities can only be characterized by their defining property as 
an image of a set of X-fermion wave functions, this minimization is equivalent to the 
original many-body eigenvalue problem on Fock space. Another problem is that the 
existence theorem's proof does not provide a construction of the energy functional. 

These two problems are circumvented by the constrained search formulation of 



Levy [O, |T3|. Consider the image of the moment map Cn = M{V{T-C)), for which 
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Q* D Cn D C. A density p is said to be A^-representable if p G Cn- Define the energy 
functional 

E : C]\f — > R, 
<t>eM i(p) 

and also Ex[p] = E[p] + {p,X) for p G Cn- 

Theorem 2. The ground state energy Eq satisfies the minimization condition 
Eq = inf Ex[p], 

P&Cn 

and the ground state density po attains this minimized energy. If p is X-representable, 
thenEx[p] = Fx(p). 

Proof: The ground state energy minimizes the expectation of the total Hamiltonian: 



Eq = min Fx{(f)) 

<t>&'H 



mm 



min (F(0) + (M(0),X)) 



= min(E[p] + (p,X)). 

pGCjv 

If </)o G 5 is a nondegenerate ground state, po G C its density, and Eq = Ex[pq] the 
ground state energy, then 

^0 = Fx{(po) = Fxipo)- 

Thus, according to the generalized Levy's theorem, it is sufficient to minimize 
the energy functional on the space of all A^-representable densities instead of all X- 
representable densities. 

To determine the universal energy functional E[p], the Lagrange multiplier theorem, 
suitably extended to Banach spaces, may be applied to eliminate the constraint to 
(f) G M~^(p) from the functional's defining equation (|^. This is achieved by regarding 
the pairing (p, X) as a Lagrange multiplier term while the operator X is the Lagrange 
multiplier itself. Suppose the vector space of operators is self-dual, i.e., the dual space 
to 0* is Q, and it is a Banach space, i.e., a norm ||X|| is defined on g. Both conditions 
are satisfied, for example, if the space is finite-dimensional. Suppose that M~^(p) is a 
submanifold of P(7i). 

Theorem 3. The following are equivalent conditions on 0o ^ M^^(p): 

(i) 00 is a critical point of restricted to M~^(p); and 

(ii) there is a G g such that (Xq, ^o) is a critical point of Fx{(p) on the space q xViH). 

The universal energy functional is given by E[p] = -F(0o) when the critical point 
corresponds to an absolute minimum. To test for a minimum, let a be a real positive 
constant, (Xo,0o) a critical point of Fx(0), and set 

F^{<P) = Fx{<P) + a\\X - X4\ (8) 
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that also has a critical point at (Xo,0o)- If the restriction of F{(f)) to M~^{p) has a 
nondegenerate minimum at (po, then for sufficiently large positive a, F^{(f)) is minimized 
at (Xo,0o)- This observation may be used to identify minima using the unconstrained 
second derivative test. The Newton method, in its Banach space generalization, may 



be applied to compute critical points for -Fx(</') See p. 219 of ^ for the theorem's 
proof. 

3. COADJOINT ORBITS 

Suppose the vector space of hermitian operators g is a Lie algebra and let G denote 
the corresponding simply connected and connected Lie group of unitary operators. The 
group action of G on Fock space naturally induces the coadjoint action Ad* of G on the 
dual space of densities. To show this, suppose (j) E H and g ■ (p E H is the result of the 
unitary representation of g E G acting on (p. If p = M^cj)), then M{g-(p) = Ad*g{p), since 

{Mig ■<!)), X) = {g-<P\Xg-<P)/{g-<P\g-^) 
= {<P\g~'Xg-<P)/{<P\<P) 
= {p,Adg-.{X)) 

= {Ad;{p),X) (9) 

for all X G 0. Thus the coadjoint action leaves the space of X-representable densities 
invariant, Ad*g : Cat — Cat. 

If G is a symmetry group of the Hamiltonian H, then the energy functional 
transforms according to the adjoint transformation, 

ExoAd* = EAd.(x). (10) 



9 



for X G and g E G. If p G C is an absolute nondegenerate minimum for Ex, 
then Ad*g{p) is an absolute nondegenerate minimum for EAd Therefore, if G is 

a symmetry group for the Hamiltonian, the coadjoint G-action leaves the space of X- 
representable densities C invariant. 

In conventional density functional theory, the abelian Lie algebra integrates to an 
abelian group of unitary operators 

G=i^g = exp (^z /(fo) , / e gj , (11) 

and G is the isotropy subgroup at each density, Ad*{p) = p. If (p E H and 
M(0) = p E Cn, then the energy functional satisfies 

E[p]<MFig-<p). (12) 

This equation provides an upper bound on the exact energy functional and may be used 
to approximate it. 

For a nonabelian Lie algebra g, the Lie group G is not typically contained in 
the isotropy subgroup Gp aX p E C^. The coadjoint orbit Op consists of all points 
Ad*Jp) E Cn as g ranges over the group G. Each orbit space Op is diffeomorphic to 
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the homogeneous coset space G/Gp. Hence the space of X-representable densities is a 
union of disjoint coadjoint orbits each diffeomorphic to a homogeneous space. 

Every Lie algebra element X defines a tangent vector to the coadjoint orbit 
surface through the iV-representable density p. This vector is the tangent to the curve 
Ixit) = Ad*g{p) where g = exp{itX), see Figure. Note that the elements of the isotropy 
subalgebra Qp correspond to null tangent vectors, and the tangent space is isomorphic to 
the coset vector space g/flp. For any two tangent vectors X,Y G g/^p to the coadjoint 
orbit through p, the bilinear form 

cUpiX,Y)=t{p,[X,Y]), (13) 

is antisymmetric, closed, and nondegenerate. Thus each coadjoint orbit is a symplectic 
manifold ||I5|. Moreover the group G acts on each orbit as a transitive group of canonical 
transformations, 

UAd;i,)iAd,iX),Ad,iY)) = ^p(X,y), (14) 
where X,Y are tangent vectors at p and Adg{X) and Adg{Y) are tangent vectors at 

Ad;{p). 

The model group theoretic Hamiltonian h may be defined now. Choose a coadjoint 
orbit Op in the dual space and let S denote the restriction of the energy functional to 
this orbit. If X is a tangent vector at p, denote the derivative of the energy in the X 
direction by 

dSp{X) = ^^£{Ad;ip))\t=o, where g = exp(ztX). (15) 

The model Hamiltonian /i G g/sp is determined by 

dSp{X) = Up{h, X) for all X G g. (16) 

There is a unique solution to this equation for /i G g/sp because the symplectic form u 
is nondegenerate on the tangent space. Note that the model Hamiltonian h depends on 
the point p in the coadjoint orbit. 

In the special case when the Hamiltonian H is an element of the Lie algebra g, the 
energy functional is just E[p] = {p,H) and the model Hamiltonian h is equal to H. 
Otherwise the model Hamiltonian h is the optimal approximation to H relative to the 
Lie algebra g in the sense that {(j)\[H — h, X] 0) = for all X G g where M(0) = p E C. 

If g is the algebra of one-body operators and the coadjoint orbit consists of 
the determinental densities, then the energy S is inherited uniquely from Fock space 
because of the one-to-one correspondence between determinental wave functions and 
determinental densities. The model Hamiltonian h on the manifold of determinental 



densities is identical to the Hartree-Fock mean field [16]. But the mean field Hamiltonian 
h is also defined here for all orbits, including those for the nondetermental densities. 
The required energy on the general orbits is the universal energy functional E[p]. 
Hartree-Fock-Bogoliubov is extended similarly by this construction using the algebra 
of quadratic quasiparticle operators and considering the energy functional on generic 
orbits of quasiparticle densities |1T7 |. 



Letter to the Editor 



8 



4. CONCLUSION 

A useful theoretical description of a many-body system is based on a wise choice for the 
operator space q. If the degrees of freedom that are most relevant to the dynamics are not 
represented in g, then the energy functional E[p] on the dual space will be intractable. 
The external Coulomb potential created by the positive charges of the massive nuclei 
dominates the physics of atomic, molecular, and solid state systems. The relatively 
weak interaction among the electrons can be treated in perturbation theory. Thus 
ordinary density functional theory, for which g is the algebra of one-body multiplication 
operators, is the appropriate theory for these interacting fermion systems. Since the 
algebra of multiplication operators is abelian, the coadjoint orbits are just points and 
no mean field approximation may be constructed in this case. 

There is no external potential acting on an atomic nucleus, and the dynamics 
is determined by a competition among short range pairing, long-range quadrupole- 
quadrupole, spin-orbit, and other hadronic internucleon forces. The relative influence of 
the different forces varies with the isotope. Hence there is no obvious choice for g that 
applies to all nuclei. If single-particle degrees of freedom are important, then the space 
of operators should be the unitary algebra u(7i) of all one-body fermion operators, i.e., 
the span of the operators b^^jba that destroy a fermion in state a and replace it with 
one in state f3. If geometrical collective modes are active, then the symplectic algebra 
sp{3, R) or one of its subalgebras is suggested O]. Since the unitary algebra u(H) and 



the symplectic algebra sp{3, R) are not abelian, a mean field approximation may be 
applied to these Lie algebras g. 
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